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ABSTRACT 

Recent high precision proper motions from the Hubble Space Telescope (HST) suggest that the 
Large and Small Magellanic Clouds (LMC and SMC, respectively) are either on their first passage 
or on an eccentric long period (>6 Gyr) orbit about the Milky Way (MW). This differs markedly 
from the canonical picture in which the Clouds travel on a quasi-periodic orbit about the MW 
(period of ~2 Gyr). Without a short period orbit about the MW, the origin of the Magellanic 
Stream, a young (1-2 Gyr old) coherent stream of HI gas that trails the Clouds ^150° across the 
sky, can no longer be attributed to stripping by MW tides and/or ram pressure stripping by MW 
halo gas. We propose an alternative formation mechanism in which material is removed by LMC 
tides acting on the SMC before the system is accreted by the MW. We demonstrate the feasibility 
and generality of this scenario using an N-body/SPH simulation with cosmologically motivated 
initial conditions constrained by the observations. Under these conditions we demonstrate that 
it is possible to explain the origin of the Magellanic Stream in a first infall scenario. This picture 
is generically applicable to any gas-rich dwarf galaxy pair infalling towards a massive host or 
interacting in isolation. 



Subject headings: 
and dynamics 
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The Magellanic Clouds are the closest known 
interacting pair of galaxies. Optical and infrared 
surveys of the system present the Clouds as two 
distinct objects separated in space by a projected 
distance of ~20 kpc. The HI distribution, how- 
ever, paints a different picture. The Clouds are 
connected by a low metallicity bridge of gas, re- 
ferred to as the Magellanic Bridge, and share a 



that the y have had a close en c ounte r in the recent 

demonstrate 



iToomre fc Toomre 



past. 

that two isolated galaxies are capable of removing 
substantial amounts of material via tides, form- 
ing pronounced features such as bridges and tails. 
Interestingly, the system also possesses a substan- 
tial trailing HI c omponent, known as th e Mag- 



ellanic Stream ( Mathewson et al. I 119741 ) . The 



Strea m is a filamentary feature of H I gas (no 



common gaseous en velope (jPutman et al. 112003 ; 



Br tins et al. 20051) . The existence of such fea- 



stars; iGuhathakurta fc Reitzel 1 119981 ) that trails 



tures suggests that the Clouds are a binary in- 
teracting pair. The Bridge in particular indicates 
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behind the Clouds f or at le ast 150° across t he sk y 
(|Braun fc Thilker I l2004t iNidever et all l201tl) . 
The Stream has historically been explained as the 
product of a tidal and/or hydrodynamic interac- 
tion between the Clouds and the MW (e.g., GN96, 



Connors et al~l 120061 : iMastropietro et al. I Eiooih 



This picture stems from the belief that the Clouds 
have traveled in an orbit that afforded multiple 
close passages between the Clouds and our Galaxy. 
However, recent HST proper motion (PM) mea- 



sureme nts of the Clouds by iKallivavalil et al. 



(|2006al |bl) (hereafter Kl and K2), independently 
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confirmed by IPiatek et al~l (|2008l ). have chal- 
lenged this picture. These studies suggest instead 
that the Clouds have, at best, completed one or- 
bit about the MW or may even be on their first 



of hig h velocity clouds by 



passage (Besla et al. 1120071 hereafter B07). This 



calls for a revised interpretation of the origin of 
the Stream. 

Kl determined a velocity of 378 ± 18 km/s 
for the LMC, which is 80 km/s higher than that 
derived from theoretical models of the Stream 
(GN96). B07 showed that a backward i ntegration 
scheme (e.g.. iMurai fe Fuiimoto II 1980f ) using the 
new velocities and an isothermal sphere model for 
the MW yields an orbit for the LMC with an apoc- 
enter >200 kpc. It is, however, unlikely that the 
rotation curve of the MW remains flat out to such 
distances (e.g., IXue et al~ 2008 ). If instead a 



more cosmologically motiva ted profile is employed 



(e.g. NFW or H ernquist, iNavarro et al. I Il996t 



Hernquist 1990a) then orbital solutions with mul- 
tiple pericentric passages are ruled out. In partic- 
ular, if an NFW model is adopted with a mass of 
1O 12 M , B07 conclude that the Clouds have just 
experienced their first close passage past the MW. 

Given the uncertainties in the adopted MW 
model, it is also possible that the Clouds have 
completed at most one orbit about our Galaxy 
within a Hubble time. This occurs if either the 
MW's mass were higher (~2x 1O 12 M ; B07), or 
if the velocity inferred from the PM measurements 
were substantially lower, e.g. if the velocity at the 
Solar circle were higher ( Shattow fc Loeb 20091: 
Reid et a"l~ll2009l ). In this case, the only previous 
pericentric passage about the Galaxy would be ^6 
Gyr ago and the apocenter of the orbit would be 
~400 kpc (i.e. larger than the virial radius of the 
MW). These values are lower limits since these 
studies assume the MW's mass is constant over 
time, whereas in the current LCDM paradigm our 
Galaxy is beli eved to have been ha lf as massive 
~8 Gyr ago (jWechsler et al. 1 120021 ). Depending 



on the mass evolution of the MW it may be im- 
possible for the Clouds to have completed multiple 
passages. 

Independent of which of the two orbital sce- 
narios outlined above is correct, there is no or- 
bital solution that brings the Clouds near the 
MW over the past 3 Gyr. However, there is 
strong evidence that the Stream is a young fea- 
ture (1-2 Gyr). Estimates of the survivability 



(2009J); iKeres fc Hernquist 



Heit sch fc Putman 
(2009) make it im- 



probable for the Stream to have survived much 
longer. The Stream also exhi bits surprisingly high 



Ha e mission ( ~750 mR; IWeiner fc Williams 



1996), which implies that cloudlets within the 



Stream are being ablated away on a 100-200 Myr 
timescale and must be contin uously replenished 
(jBland-Hawthorn et al. 1 120071 ). The lifetime of 
the Stream poses a problem for all past numerical 
models, which invoke some combination of MW 
tides and/or ram pressure stripping to form the 
Stream. These models require at least one com- 
plete orbit about the MW, implying an incompat- 
ible age of at least 6 Gyr. 

As such, regardless of whether the Clouds are 
on their first or second passage about the MW, 
in the context of the origin of the Stream we are 
left with the same problem: How can the Stream 
have formed without a complete orbit about the 
MW? Or, more generically, how can a pronounced 
tail be formed from a pair of dwarfs on their first 
infall towards a massive host? 

Based on the above considerations, we explore 
the following scenario using simulations. We as- 
sume that the Clouds were a stable binary sys- 
tem and adopt the simplifying hypothesis that the 
Clouds are on their first passage about the MW in 
order to illustrate how the Stream can form with- 
out relying on a close encounter with the MW. We 
postulate that the Magellanic Stream and Bridge 
are in fact a classical bridge and tail caused by the 
tidal interaction between the Clouds before they 
have been accreted by the MW. The MW poten- 
tial shapes the orbit of the Clouds and thereby 
controls the appearance of the tail, causing the 
line-of-sight velocities and spatial location of the 
tail to be as observed in the Stream today. 

Given these assumptions we attempt to explain 
the following observed features of the Stream: 1) 
The absence of stars; 2) The ~150° extent; 3) The 
spatial location projected on the plane of the sky; 
4) The line-of-sight velocities along its length; 5) 
The HI column densities; and 6) The pronounced 
asymmetry between the trailing and leading com- 
ponents. 
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2. Methodology 

Simulations were carried out using the N- 
body/SPH code GADGET2 (Springel et al 2005). 
Star formation is not included for simplicity, but 
the gas is allowed to cool radiatively. 

We model the MW as a static NFW potential 
of mass of 1.5 x 10 12 M Q , i? vir =240 kpc and a 
concentration parameter of 12. Dynamical friction 
from the MW halo is not explicitly accounted for, 
but is expected to have little impact on the orbit 
in a first passage (see B07, figure 4). The L/SMC 
are both modeled using Hernquist profiles for their 
dark matter content, and exponential gaseous and 
stellar disks, where the scale length of the gas disk 
is 6 (LMC) or 5 (SMC) times that of the stellar 
disk. Such extended gaseou s disks are not atypica l 
for isolated dwarf galaxies (|Swaters et al. 1120021 ) . 



lisions between the Clouds have caused material 
to be lo osely bound to the Clo uds is not novel 



The stellar and gas masses for the Clouds are 
well constrained within their respective observ- 
able limits. However, the total dark matter con- 
tent of these galaxies is unknown. All previous 
models of the Magellanic system have assumed 
the L MC is tidally truncated to a radius of 15 



kpc (jvan der Marel et al. 1 120021 ) , resulting in a 
total mass estimate of 2-3 xlO 10 ^ (e.g., GN96 , 

Bekki fc Chibal [2005). 



Murai fc Fuiimoto I fl98cl 



On an orbit where MW tides are largely incon- 
sequential, the L/SMC will not be truncated. In- 
stead, we use current halo occupation models to 
relate the observed stellar mass of the LMC to its 
original halo mass before infall into the MW halo 
(Ivan den Bosch et al. 1 l2003t IConrcw et al. 



2006; 



Wang et al. 20061 ). In Table Q]we summarize our 



adopted model properties. The L/SMC are found 
to have infall masses an order of magnitude larger 
than employed in previous models. 

3. Simulation Results and Comparisons 
with Data 

We believe that the key to understanding the 
origin of the Stream lies in understanding the in- 
teraction history of the Clouds themselves. The 
majority of previous models of the Stream have 
assumed that the SMC is in a circular orbit about 
the LMC (separation of ^20 kpc). Using the cos- 
mologically expected infall masses, the dynamical 
friction timescale for such an orbit would be much 
less than a Hubble time, so the SMC's orbit about 
the LMC cannot be circular. The idea that col- 



(GN96, Hel ler fc Rohlfs 1 119941 ) . However, the 



morphology resulting from collisions between the 
Clouds as the SMC travels in a highly eccentric or- 
bit about the LMC, independent of the MW, has 
never been explored. 

We have used the current projected separation 
between the Clouds ( 23 kpc), and the observed 
relative velocity (~ 105 ± 42 km/s; K2) to con- 
strain the orbital history of the SMC about the 
LMC. The resulting orbit is nearly parabolic (e= 
0.7), with an apocenter of ~ 100 kpc (see Figure 
Q]). The SMC disk is oriented 90° with respect to 
its orbital plane about the LMC. If the SMC were 
in a coplanar, retrograde orbit about the LMC, no 
material would be removed with this orbital con- 
figuration. This implies that dispersion supported 
material (e.g. dark matter or stellar halo) will be 
unaffected within the SMC's disk radius. 

The presented orbital solution is not unique; 
however, a highly eccentric orbit is required to pre- 
vent the Clouds from merging. Moreover, these or- 



bits are cosmologic ally typical (see, e.g.. iBenson 



20051 : IWetzel IlioTol) . Any eccentric orbit that al- 



lows for high speed encounters between the Clouds 
will yield similar bridge and tail structures. 

We simulated the interaction between the 
Clouds as the SMC travels along the orbit shown 
in Figure 1 (top left) starting > 6 Gyr ago. At 
each close passage between the Clouds, gas is re- 
moved from the SMC by LMC tides, forming a 
tidal tail and bridge. However, as the Clouds 
move apart, the tidal material falls back towards 
the L/SMC disks. These features are thus tran- 
sient phenomena. Little material is removed from 
the LMC, despite its extended gaseous disk com- 
ponent: SMC tides are ineffective. 

The pericentric passage at 4 Gyr in Figure Q] 
(top left) results in significant gas removal from 
the outer regions of the SMC's extended gas disk, 
whereas the smaller stellar disk remains intact. As 
the SMC travels away from the LMC, material 
that was tidally removed by the LMC stretches 
our to distances as large as 100 kpc from the 
Clouds. This material will eventually form the 
Stream (Figure [TJ bottom right). Hence, an ex- 
tended gaseous tail is produced without the aid of 
MW tides or ram pressure. 
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We assume that the Clouds have been an inter- 
acting pair for a significant fraction of a Hubble 
time and have crossed the virial radius (240 kpc) 
of our Galaxy ~1 Gyr ago. From Figure [TJ a 100 
kpc long tail is formed after ~5 Gyr. We there- 
fore stop the simulation 1 Gyr earlier (i.e. at the 
timestep corresponding to the bottom left panel 
of Figure [T]) and place the binary system outside 
the virial radius of the MW. We then allow the 
Clouds to travel to their current observed loca- 
tions on an orbit consistent within la of the PMs 



oftheLMC (Kl. lPiateket al. Il2008h . We did not 



attempt to reproduce the SMC PM determined by 



K2 and iPiatek et al. I ([20081 ) since they are dis- 
crepant. However, the resulting SMC line-of-sight 
velocity agrees well with the observed value. 

In Figure [5] we show the resulting stellar dis- 
tribution (top) and HI gas column density map 
(middle) for our simulated stream: there are no 
observable stars in the simulated stream and the 
stream extends ^150°. The white line indicates 
the current location of the observed Stream and is 
well matched by the simulation. Notice that the 
past orbits (yellow lines) are not co-located with 
the simulated stream, as expected since the north 
component of the LMC PM vector is not aligned 
with the Stream (see B07, figure 9). This spatial 
mismatch is a natural result of our model: mate- 
rial is removed from the SMC along the LMC-SMC 
binary orbital plane, which is not co-planar with 
the orbital plane of the Clouds about the MW. 
The SMC disk was oriented 90° to the L/SMC bi- 
nary plane in order to maximize this offset. This 
further implies that the SMC is seen ~edge on 
from our viewing perspective, explaining its sur- 
prisi ngly large observed line-of-sight dept h (~5 
kpc; ISubramanian fe Subramanian 20091) . 

Figure [2] (bottom) illustrates the resulting line- 
of-sight velocities for the simulat ed stream. The 
white line is a fit to the data ( Putman et al. 



2003: Nidever et al. 



2003) [hereafter, N08] and 
shows good agreement with the simulation. The 
line-of-sight velocities along the new orbits (yellow 
lines) are much larger than those observed along 
the Stream: since the Stream and the orbits are 
not co-located, their line-of-sight velocities are not 
similar (see also B07, figure 20). The MWs gravi- 
tational field serves to stretch the stream and mod- 
ify its galactocentric distance: the tip of the simu- 
lated stream is located ^100-140 kpc away. Using 



these distance estimates, we find that the mass in 
our simulated stream matches th e observed value 
(4.8 x lO 8 M : lBriins et al. 1^0051 ) to within a fac- 
tor of 2. 

The simulated gas column densities range from 
10 18 — 10 21 /cm 2 as observed, although the exact 
column density gradient along the length of the 
Stream is not reproduced. The column density 
is not homogeneous across the width of the simu- 
lated stream: the inclusion of metal cooling, ion- 
ization and confinement/interaction by /with the 
ambient MW halo gas will likely aid in reproduc- 
ing the bifurcated, filamentary nature of the ob- 
served Stream. 

The leading component of the simulated stream 
is much smaller than the trailing component, as 
observed. This arises because the leading tidal 
arm from the SMC falls towards the LMC, while 
the tail continues to grow (Figure Q]). However, 
the leading component in the model is not in the 
correct location on the sky when compared to ob- 
servations. This material leads the orbit, which is 
not aligned with the Stream. Consequently, this 
apparent problem w ill also occur in t he traditional 
tidal m o dels of e.g. iRuzicka et al. I (120091) :lB ckki 



2008 ): IConnors et al. I (2006): Bek ki fe Chiba 



([20051) and iGardiner fc Noguchil (|l996l hereafter. 
GN96). This indicates that hydrodynamic pro- 



cesse s such as ram pressure (Ma stropietro et al 



2005) are needed to shape the final appearance of 
the simulated stream. 

4. Discussion and Conclusions 

We have shown that tidal interactions between 
the Clouds are sufficient to remove a substantial 
amount of material from the SMC without the 
aid of MW tides or hydrodynamic interactions. 
Our goal was not to reproduce every detail of 
the Stream, but merely to show that a plausi- 
ble Stream model can be formed without a pre- 
vious passage about the MW. Nevertheless, many 
features of the simulation fit the data remarkably 
well. 

We have explained the absence of an observ- 
able (brighter than Vband 27 mag/arcsec 2 ) stellar 
counterpart to the Stream by placing the SMC on 
an orbit about the LMC with a large impact pa- 
rameter, which allows its compact stellar disk to 
remain intact (see also, lYoshizawa fc Noguchi 
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20031 ). The 150° extent of the Stream is repro- 
duced owing to the eccentric orbit of the SMC 
about the LMC (apocenter -100-50 kpc) and MW 
tides after the system has been accreted. Such an 
eccentric orbit is likely only viable in a first in- 
fall scenario, as MW tides would have disrupted 
the system at the previous pericentric approach. 
The spatial location of and velocity gradient along 
the simulated stream are correct and offset from 
those of the orbits: material is removed in the 
L/SMC binary plane, which is not coplanar with 
the orbit of the Clouds about the MW. The simu- 
lated HI column densities also match the observa- 
tions, although the exact column density gradient 
is not reproduced without additional physics. Fi- 
nally, the strong asymmetry between the leading 
and trailing stream components is a natural conse- 
quence of a classical tidal bridge and tail scenario. 

There are, however, other features that require 
additional physics to address, such as the claim of 
correlated bursts of star fo rmation in both Clouds 
([Harris fc Zaritskv I [20091) . In our new scenario, 



correlated bursts of star formation may correspond 
to close passages between the Clouds rather than 
with the MW. Furthermore, without ram pres- 
sure and stellar feedback we cannot explore the 
claim by N08 that half of the Stream originates 
from the LMC as a stellar outflow (see also, B07 



and lOlano 1 12004). We intend to examine these 



issues and more detailed comparisons of our simu- 
lations to the Magellanic system in future studies. 
However, we note that recent me tallicity measure- 
ments at the tip of the Stream ( Fox et aT~ll2010l ) 
indicate that the Stream is extremely metal poor 
and thus inconsist ent wi th a stellar outflow sce- 

also find that the oxy- 



Fox et al. 



nario. 

gen abundance at the tip is more consistent with 
that of the SMC rather than the LMC. As such, a 
model where the Stream originates primarily from 
the SMC is not ruled out by observations of the 
Stream. 

The simulation results presented here with re- 
spect to the L/SMC interaction have broader 
applicability than just to the Magellanic sys- 
tem. Within the current LCDM paradigm ha- 
los at all scales are expected to build up their 
mass hierarchically. Interacting dwarf galaxies are 
therefore cosmologically expected both in isola- 
tion and within larger halos ( Simha et al. 2010t 
D'Onghia et al. 1 120091 : [Knebe et al. II2006I ). The 



isolated interacting Magellanic type galaxies NGC 
4490/85, w hich are surrounded by an extended HI 
envelope ( Clemens et al. 19981) . and the inter - 
acting M51/NGC 5195 pair dHernauist I Il990bl) 
may be observational analogs of our initial L/SMC 
system. The presented model thus illustrates that 
dwarf-dwarf galaxy tidal interactions are a pow- 
erful mechanism to morphologically change dwarf 
galaxies without the need for repeated interactions 
with a massive host. 
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Fig. 1 — The orbit of the SMC about the LMC (without the MW). Top left: The separation between the 
Clouds as a function of time. The chosen orbit is highly eccentric (e=0.7) and decays rapidly owing to 
dynamical friction. The Stream forms at 4 Gyr (red box) and the time today would correspond to ^5.2 Gyr 
(blue box). Subsequent Panels: The HI gas column density is plotted as contours over the stellar distribution 
for specific points in the orbit. Gas contours span a range of 10 18 — 10 20 /cm 2 , where each contour represents 
an increase in column density by a factor of 1.5. HI gas is identified as gas at temperatures below 12000 K 
and column densities larger than 10 18 /cm 2 , although a background ionizing field is not included. A bridge of 
gas connecting the Clouds and a 100 kpc long gaseous stream with no stellar counterpart is formed without 
the aid of MW tides or ram pressure. 
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Magellanic Longitude 



Fig. 2. — The stellar surface brightness, HI gas column densities and line-of-sight velocities of the simulated 
Magellanic system. Top panel: The resulting stellar distribution is projected in Magellanic coordinates 
(N08), a variation of the Galactic coordinate system where the Stream is straight. The distribution is color- 
coded in terms of Vband surface brightness. The past orbit of the L/SMC are indicated by the blue lines. 
Middle panel: The HI gas column densities of the simulated stream range from 10 18 — 10 21 /cm 2 , as expected 



( Putman et al. 1120031 ). The white circle indicates the observed extent of the LMC's HI disk: the simulated 



LMC is more extended than observed, indicating ram pressure likely plays a role to truncate the disk. In 
both the top and middle panels, the solid white line indicates the past orbit of the SMC according to the 
old theoretically derived PMs (GN96) which was a priori chosen to trace the Stream on the plane of the 
sky. The true orbits (determined by all PM measurements) for the L/SMC are indicated by the yellow lines. 
Bottom panel: The line-of-sight velocities along the simulated stream are plotted and color-coded based on 
HI column density, as in the middle panel. The white line is a fit to the observed data (N08). The LMC 
disk is too extended, causing a larger velocity spread than observed. The line-of-sight velocities along the 
past orbits of the L/SMC are indicated by the yellow lines, which do not follow the true velocities along the 
Stream (e.g. B07 figure 20). The Stream is kinematically distinct from the orbits of the Clouds. 
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Table 1 

L/SMC Initial Conditions Adopted 



Property LMC SMC 



M* (M Q ) 


2.2 x 10 9 


1.3 x 10 8 


M gas (M ) 


1.4 x 10 9 


1.1 x 10 9 


M halo (M ) a 


1.8 x 10 n M© 


2.5 x 10 10 


V200 (km/s) 


82 


42 


Concentration 


9 


15 


Stellar scale length (kpc) 


1.7 


0.7 


Gas scale length (kpc) 


10.2 


3.6 


Nstars 


100000 


100000 


Ngas 


300000 


300000 


Nhalo 


100000 


100000 



a The total halo mass is determined using the observed 
stellar m ass of the LMC(SMC) M* = 3 x 10 9 M W (3 x 



10 s M y) (|van der Marel et al _ 
2004h and the relations from [Wang et ah 



2002t Stanimirovic et al. 



(120061k Note 
is less than that quoted in this table 



that the observed M 
in order to account for star formation in future studies 
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